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Abstract. Multifractal measures are applied to three IRAS galaxy subsamples selected by colour from the PSCz catalogue. As 
shown by a generalised dimension spectrum, hot IRAS galaxies are found to be less clustered than cold galaxies, but the differ- 
ence is very small. An alternative tool, the conditional multifractal dimension spectrum reveals an apparent peculiarity of the 
distribution of hot galaxies, especially at high orders. The conditional multifractal measure basically measures the environment 
of selected galaxies. A detailed analysis of the distribution of galaxies with their number of neighbours shows that cold galaxies 
are more likely to be in over-dense regions than hot galaxies. Further studies suggest that since the colour of IRAS galaxies is 
a good trace of star formation rate, it is possible that we have statistical evidence here for enhanced star formation rate due to 
galaxy interactions. 
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1. Introduction 

Galaxies are inevitably affected by their environment. It has 
been found that early-type galaxies are far more common than 
spirals in centres of rich clusters, while the opposit e is true in 
other parts of the universe (.Dressier. .1980: .Caon & Einastol 
Il995t IPressler et all 119971) . There is the same phenomenon 
in galaxies' spatial distribution: early-type ga laxies ar e more 
strongly clustered th an late-tvp es ^e.g .Hermit et all Il996t 
IShepherd et aP. 1200 it iMadgwic k et all l2003b . and high lumi- 
nosity galaxi es are more strongly clustere d than lo w luminos- 
ity oiies (e.g..lHermit et al.[ll996tlLin et al.L.1996:.Guzzo et all 



3. 



l2000t iNorberg et all 1200 ih . It is clear that galaxies are biased 
tracers of the mass, which has invoked the concept of bias to 
make connections between statistics of galaxies and the mass 
(e.g.. iKaiseil 1 1984 iBardeen et all Il986[ IMo & Whit^ Il996t 
IPekel & Laha v. 1999). Galaxy samples selected by different 
criter ia are frequently studied to quantify t he phenomenon 
(e.g.. lNorberg et alll2002llZehavi et alll2005l) . 

Among the many galaxy surveys, the IRAS galaxy sample 
plays an important role in studying cosmic large-scale struc- 
ture and galaxies' formation and evolution. Observationally, 
it affords large sky coverage, uniform flux calibration, good 
position accuracy, and insignificant galactic absorption. IRAS 
galaxies usually are good objects of strong star formation activ- 
ities. The majority of galaxies observed in the infrared band are 
spirals, which tend to avoid rich galaxy clusters, so in general 
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the infrared selected galaxy samples have a lower clustering 
amplitude. 

The infrared emission from galaxies are produced by their 
gas (dust) component, illuminate d main l y by UV emission 
from stars. iRow an-Robinson & Cr awford! (119891) decomposed 
the spectrum of IRAS galaxies into three components: (1) a 
cool disk component emission from interstellar dust activated 
by galaxy starlight; (2) a warmer "starburst" component in- 
terpreted as a burst of star formation in the galaxy nucleus, 
whose spectrum is well fitted by a model of hot stars embedded 
in optically thick dust cloud; and (3) a "Seyfert" component 
that originates in a power-law continuum source within a dust 
cloud related to the narrow-line region of the compact source. 
The peak of the disk component's radiation is at the 100//m 
and the "starburst" part reaches maximum at 60yL(m, while the 
"Seyfert" component's peaks are at 12/im and 25jum; thus, the 
far-infrared radiation from the IRAS galaxies is a composition 
of ambient cirrus emission and localised emission from regions 
of active star formation. The relative contribution of the two 
components may be quantified by the dust temperature inferred 
from the flux ratio S 100/5 60- 

Therefore, SimJSa), or equivalently the dust tempera- 
ture, is a good measure of the star formation rate (SFR ) 
of IRAS galaxies JSaunders et all Il990l: iMann et all Il996l) . 
There are arguments that in deep gravity potential wells, S FRs 
of galaxies are effectively suppressed (e.g . , Young et al.. 1996t 
iBalogh et all Il998l: iBIanton et all Il999l) ; nafiirally, the segre- 
gation phenomenon in clustering is expected in the temperature 
or the SFR selected IRAS galaxy subsamples, which is sup- 
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ported by the observation that late-type optical galaxies with 
higher SFRs have smaller clustering amplitudes than early-type 
galaxies. ^^^_^^^ 

SurprisinglvJ Mann et alJ (119961) found an unusual cluster- 
ing scenario of their "warm" and "cold" subsamples of QDOT 
divided by the 36K temperature criteria. Their "warm" galax- 
ies are more strongly clustered than the "cool" galaxies, which 
is oppo site to what one would expect. However, Hawkins et alJ 
11200 11) split the PSCz catalogue, which is superior to QDOT, 
into three subsamples, "hot", "warm", and "cold", which fall 
into temperature ranges > 3 IK, 28K~ 3 IK, and < 28K, re- 
spectively, and they found that the hot galaxies are actually 
clustered less strongly than those cold galaxies on scales be- 
tween 1/i^'Mpc and 10/;"'Mpc, though the trend is weak. 

Here, the problem of the colour dependence of the IRAS 
galaxies' clustering is revisited using a multifractal tool. For 
many years, fractal dimensions were used to test if the distri- 
bution of galaxies in the universe is homogeneous. Actually, 
it could have more applications in modern statistical cosmol- 
ogy than simply as a challenge to the Cosmological Principle. 
The benefit of using multifractal is that we have a whole spec- 
trum of dimensions to describe the distribution. Actually, in 
the past, t he fractal was already used to study morphology seg- 
rega tion (Wenetal., 1989; Dominguez-Tenreiro & Martinez, 
11989.: .Dominguez-Tenreiro et aL .1994: .Best et aL .1996). In 
addition to conventional fractal dimensions, a new tool, con- 
ditional dimension, was designed to measure differences be- 
tween two samples. The idea of conditional dimensions was 
inspired by concepts of conditional and relative multifractal 
spectra (Rich ards & Scheuringlll997l) . Details of the subsam- 
ples' construction are given in Sec. 2. In Sec. 3 we present re- 
sults from ordinary multifractal measurement; Sec. 4 is ded- 
icated to new statistics, the conditional multifractal measure. 
Conclusions and discussions are in the last section. 



2. Sample construction 

2.1. The main sample from PSCz 

The PSCz catalogue contains 15411 ga laxies covering 84% of 
the sky with 60yum flux, Sm > 0.60 dSaunders et all I2OOO); 
redshifts are available for 14677 galaxies with a redshift me- 
dian of ~ 0.028. Those galaxies with galactic latitude ^ < 10° 
are thrown away to exclude the high galactic absorption plane. 
Further selection criteria are applied by the galaxy's co-moving 
distance calculated with Ho = IOO/1 'kms ' and Q = 1, which 
is 10/1 'Mpc < r < 250/z"'Mpc to avoid the Galaxy's local in- 
fluence and make the sample in the remote end not too sparse. 
The final main sample for this work has 1 1463 galaxies. 

We find that the colour of galaxies shows an apparent radial 
gradient. It is known that IRAS galaxies experience rapid lumi- 
nosity evolution with redshift; we need to che ck if the colour 
radial gradient is an evolution efi'ect. At 60um. ISaunders et alJ 
( 1199 0) gives Leo oc (1 h- 2)^=^'. Although the exact luminos- 
ity evolution at 100/im is not clear, we can use the luminosity 
evolu tion at 90;um, which is Lgo oc (1 -h z)^ '*^' dSerieant t 
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120041) . as an approximation, so the colour is roughly propor- 
tional to (1 + z)°"^. The largest redshift of the sample is 0.089 



Fig. 1. Redshift distributions of galaxies of the three subsam- 
ples. 



which introduces correction to the colour by a factor of ~ 1 .035 
for the farthest galaxies. Even by this factor, since the mean 
colour of those remote galaxies is ~ 1 . 1 of our sample, the 
modification to the colour is about 0.04, which is relatively 
small. In fact, after applying such correction to the colour, 
changes to the straight line fitting the \og{S \oo/S(^)-\og r rela- 
tion are negligible. The colour undergoes very weak evolution 
in the redshift regime of the sample in analysis here. 

Because the main sample is flux-limited, those faint galax- 
ies do not enter the list, and the mean luminosity will obviously 
increase with distance. Actually, there is a correlation between 
luminosity and colour: less luminous galaxies are likely to have 
higher colours. Therefore, those galaxies with high colour are 
unlikely to be selected in our main sample. The mean colour 
is low in places at a large distance, while the radial grad ient 
of colo ur of the main sample is a selection efifect. Mann e t alJ 
(11996) applied a similar mean colour-distance relation to the 
construction of subsamples from QDOT because they thought 
such a radial gradient of temperature (colour) reflects the intrin- 
sic evolution. They found that after this coiTection, the cluster- 
ing diff'erence between "warm" and "cool" subsamples disap- 
pears. Now it is clear that their procedure of correction actually 
makes the subsamples quite contaminated: the nearby galaxies 
of high 5 100/5 60 are likely to be classified as "warm", while 
those remote galaxies of low colour are marked as "cool". 



2.2. Colour subsamples 

Colour subsamples a re generated using the method of 
iHawkins et al.l ( 1200 11) . "Cold" galaxies are defined as 
S\Qo/Seo ^ 2.3, the "warm" subsample consists of those of 
1.7 < S ioo/S(a < 2.3, and the hot galaxies are the remaining 
galaxies with Sm)/S(,o < 1.7. Finally, there are 3917, 4010, 
and 3536 galaxies in the "cold", "warm" and "hot" subsam- 
ples, r espectively. Differences in numbers with iHawkins et alJ 
(120011) are due to our distance cuts. The redshift distributions 
of galaxies of these colour subsamples are in Fig. [2 
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To estimate the error bars, we also created 20 bootstrap sub- 
samples for each colour subsample, as well as the main sample. 
Scatters at the 1-cr level over the 20 bootstrap subsamples are 
taken as our errors. 

3. Multifractal analysis 

3.1. The statistics and selection function 

The multifractal measure in use is constructed from the parti- 
tion function 
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with p(i) - ni(r)/N, where n,(r) is the count of objects in the 
cell of radius r, centred upon an object labelled by / (not in- 
cluded in the count), after corrections for the boundary effect 
and selection function ip{r) {ijj- 1 for volume-limited sample). 
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and fi{r) is the volume fraction of the sphere centred on the 
object of radius r within the boundaries of the sample. From 
previous works, we know that on scales larger than about 
20 ~ 30/i~'Mpc, th e distributions of galaxies are homogeneous 
( iPan & ColesL 1200 0): the comparison here focuses on scales 
less than 10/z"'Mpc to ensure reliable scaling features. Since 
the boundary correction does not significantly mo dify the result 
for th e PSCz sample, especially on small scales (IP an & ColesL 
I2OO2I) . to reduce the computation burden, the volume correc- 
tion method is adopted to deal with boundaries. 

The Renyi dimension Dg is derived from the mass exponent 
T{q) by 



Dg = T(q)/(q - 1) . 



(4) 



When q - 2, the partition function Z is the correlation in- 
tegral and D2 is the correlation dimension. For each value 
of q, Z(q) and Dq give information about the scaling proper- 
ties of different aspects of the density field. For high q, they 
are dominated by high-density regions, while for low q and 
q < Q the measure is weighted toward low-density regions. 
Note for q - \, the partition funct ion is defined differently as 
Z{q = l,r) = i;,log/:'(0 ~ r^' (cf. 'Martin ez & Saaill2002l) . In 
this paper, we concentrate on the moments of order q>2. 

From Fig. n] we know that the redshift distributions of the 
galaxies of the subsamples are not the same. Selection func- 
tions of subsamples have to be estimated separately. Of course, 
it is optimal to use a maximum-likelihood method to gener- 
ate selection functions for all subsamples. Here we just simply 
smooth those curves and then fit the distance-density relation 
via a nonl inear least-square fittin g to a double power-law func- 
tion, as in lSaunders et alJ yOOy), 

1-03/r) 
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in which i/r,, a, r,, 7, andyS are fitting parameters. A comparison 
of the partitio n functions of the mai n sample with the selection 
function from Saun ders et al.l ( 12000 ) indicates the accuracy of 
this "naive" method and that our simple fits are good enough 
for the statistics (see Fig.|3- So, for all colour subsamples, we 
adopt the nonlinear least-square fitted selection functions. 

3.2. Results of multifractal analysis 

Measurements of Z(q, r) at different orders of our colour sub- 
samples are displayed in Fig.|3j error bars are estimated from 
20 bootstrap subsamples for each colour subsample. It is obvi- 
ous that on large scales the distributions of all subsamples be- 
come homogeneous, which means the generalised dimensions 
are all close to 3. The scales of interests are therefore below 
~ 20/1" 'Mpc, a level at which the distribution is far from ho- 
mogeneity. On the other hand, to avoid severe discreteness ef- 
fects, measurements on scales of less than l/i^'Mpc are also 
cut off. In the regime of ~ 1 - 10/i^'Mpc, local dimensions 
are approximately constant with small fluctuation, and the gen- 
eralised dimensions Dq are calculated by fitting the partition 
function in this scale range. Note that Z(q, r), on small scales 
of all the colour subsamples, appears irregular with large error 
bars when ^ > 3, so one should be careful to use the dimen- 
sionality obtained by a regression of log Z{q > 3) - log r. 

Comparisons of the Dq of the subsamples are displayed in 
Fig El The most important dimension is the correlation dimen- 
sion D2. For the "hot", "warm", and "cold" subsamples, D2 is 
2.34+0.07,2. 16+0.06, and 2.14+0.03,respectively. Therefore, 
cold galaxies are distributed with a smaller fractal dimension 
than hot galaxies, which tells us that cold galaxies are more 
strongly clustered. 

The dimensionality differences among the three subsam- 
ples become smaller for higher order moments. High order 
moments are dominated by those cells from regions with large 
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Fig. 2. Dimension spectrum of ^ e [2,3] of the main sample 
under correctio ns of different selection functions. One is from 
ISaunders et alJ ( 12000 ) and the other is acquired from a direct 
numerical fit of N-z- The generalised dimensions are fitted at 
scales r < 10/i^'Mpc. 
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Fig. 3. Z{q,r) of q=2, 3, 4, and 5 for colour subsamples. Only error bars of the warm subsample are plotted; the other two 
subsamples have approximately the same size error bars. 
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Fig. 4. Generalised dimensions Dq of q e [2,3] in r E 
(1.5, 10/z^'Mpc) for colour subsamples; points of the "warm" 
and the "hot" subsamples are slightly shifted to the right for 
better presentation. 



number of galaxies, while contributions from cells in sparse re- 
gions is effectively suppressed. So, it is possible that the richest 



regions of all the subsample obey nearly the same scaling law. 
Unfortunately, the method we used here is unable to estimate 
the other half of the dimension spectrum of ^ < 2, otherwise, 
it is probable that we would find bigger differences. Of course, 
due to large uncertainties, one should be conservative about this 
observation. 

Although the partition functions of all three subsamples 
seem to asymptotically agree with each other on large scales 
within error bars, their selection functions and corresponding 
normalisation factors are different, so we cannot compare the 
amplitudes of the partition functions directly. 



4. Conditional multifractal measure 

4.1. Definition 

It is a tradition to split a main sample into several subsam- 
ples to explore differences in their spatial distributions. This 
treatment wastes information provided by the configuration 
of relative positions among galaxies of different subsamples. 
For many years the only tool to overcome this was the cross- 
correlation function. The mark correlation functions, an ele- 
gant measure recently introduced into statistical cosmology are 
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ideal for quantifying the morphology dependence of cluster- 
ing (Beisbart & Kerscher, 2000; Szapudi et al., 2000; ShetJ 
I2OO5I) . Using the mark correlation functions it has already been 
found that there is no luminosity d ependence of clustering in 
the PSCz catalogue JSzapudietaU EoOO). 

The environment of galaxies of a particular type is a con- 
cept closely related to cross-coiTelation. The partition sum of 
Eq.^measures moment of the counts of neighbours; it is natu- 
rally a good tool to study the environment, if we include galax- 
ies of all types as neighbours instead of only galaxies of the 
same type as the centre galaxy. The new measure is given 
a nam e of "conditional multifractal" ( iRichards & Scheurina 
I1997I) . In the following, we give its formal definition. 

Let the main sample of A^ points be 5 = U^^^S„„ of which 
Sm=\,...,M is a subsample in which there are N^ objects marked 
by their position vectors r„,, and{rj\j - 1,...,N} - U^^j{r,„,|/ = 
1, ...,Nm}- The conditional partition function of Sm is thus 

*P(|r,-r,„,|-r)" 



^^"''(^'^)=fZ 
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For our flux limited sample, we need to apply the same coiTec- 
tions to the kernel that we did in Eq.|2 

Formally, by writing the measure of neighbours in the /'' 
cell as ^j, i - 1, ...,A^, the conditional mass exponent r^ and 
the spectrum of the generalised conditional dimension Ddq) of 
the subsample S,,, is 
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When q - 2, Zdq - 2, r) is the conditional correlation integral 
and Dc(2) is the conditional correlation dimension. 

The subset is a collection of specially picked sampling 
points of the main set, since 
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the relation between Dg and Ddq) is 



Da 



M 

z 

m=l 



q-l 



2j ,= 1 A'm, 



y™ n' 
Lj=i f^j 



q-l 



X D'^"\q) 



IS an 



uniform dilution 



(9) 



of S, 



If a subsample S,, 

ZS'l/^mTVSjLijuf' = N,„/N and oriq) = D^. Any 
deviation from such uniformity is embedded in the conditional 
dimension; the conditional multifractal measure is therefore an 
effective tracer of differences among subsamples. 



4.2. Results of conditional multifractal analysis 

Conditional partition functions and conditional generalised di- 
mensions are calculated for each colour subsample, and results 
are displayed in Figs.|5]and|6l Recall that neighbours of an ob- 
ject are thos e in the main sample, s o only the selection function 
provided by Saunders et al. KOOOl) is needed. The error bars are 
from 20 bootstrap subsamples. 
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Fig. 6. 
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The conditional dimension spectrum within q e [2, 3] 
€ (1, 10)/!"'Mpc. Points of the "warm" and the "hot" 
subsamples are slightly shifted to the right. The solid line with 
error bars is the dimension spectrum of the main sample. 



Contrary to the ordinary multifractal measure, the condi- 
tional dimensions 0^(2) of the colour subsamples are consis- 
tent with each other within error bars, while higher orders dif- 
fer. We can see that conditional dimension spectrum of the 
"hot" and the "warm" subsamples are in good agreement with 
the spectrum of the main sample within error bars, which 
means the two subsamples are basically uniformly diluted real- 
isations of the main sample. The "cold" subsample has a quite 
different dimension spectrum; its conditional dimensions are 
larger than those of the "hot" subsample at high orders, e.g., for 
q=3, D^'^"'''* = 1.70 + 0.04, and of"'^ = 1.50 + 0.08. One thing 
that needs to be addressed is that while the conditional dimen- 
sions are a measure of the environment of the central object, a 
larger conditional dimension does not mean that the environ- 
ment is less clustered, it just tells us that the class of objects 
that are measured are in regions of a special scaling property or 
"strangeness", in fractal language. 

On large scales, the distribution of galaxies is asymptot- 
ically homogeneous, and the Zdq, r) of different subsamples 
are consistent within error bars (Fig. |5}- So, if the averaged 
environment of subsamples on small scales is different, it will 
be shown by the amplitudes of Zdq, r). For q=2, Z^, the mean 
number of neighbours, of all the colour subsamples agree with 
each other within error bars. When q > 2, zf" is much 
smaller than Z). ■ Since {^i)coid - itAhot, the variance in the 
number of neighbours around cold galaxies ((// - {ii))^)coid - 



y(cold), _ 



y{cold)2 



Z^ {q = T) — Z). (q — 2) is smaller than that around the 
hot and the warm galaxies, hence the skewness. 

At high q, Zdq, r) has similar properties to the ordinary 
Z{q, r) and is also weighted toward rich cells. It is interesting to 
see the distribution of richness around galaxies. As an example, 
the distribution function of the number of neighbours around 
galaxies of colour subsamples at scale r ~ 4/i^'Mpc are shown 
in Fig.0 Clearly, hot and warm galaxies have higher possibil- 
ities of being located in very sparse regions than cold galaxies, 
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Fig. 5. Zciq, r) of q=2, 3, 4, and 5 of the PSCz colour subsamples. Error bars of the "warm" subsample are not plotted; they are 
very similar in size to those of the "hot" subsample. 



and cold galaxies are more likely distributed in rich environ- 
ments. Cases at other scales less than 10/z^'Mpc are similar. 

5. Discussion 

In this work, colour subsamples of PSCz galaxies are analysed 
with two kinds of multifractal measures; all analysis clearly 
denotes that their spatial distributions are different. The or- 
dinary multifractal measure of our three colour subsamples 
gives results consistent with the normal morphological segre- 
gation scenario of galaxy clustering. Hot galaxies are clustered 
less strongly with a larger correlation dimension D2 than cold 
galaxies. The difference is not very strong, and actually, higher- 
order dimensions fail to detect a significant difference. 

Details of the spatial distribution of galaxies of different 
colours are more clearly displayed by the new statistics — the 
conditional multifractal that is designed to measure the envi- 
ronment around galaxies of a particular type. The conditional 
correlation integrals and conditional coiTelation dimensions are 
the same for all colour subsamples within error bars; at high or- 
ders q > 2, conditional multifractal dimensions detect marginal 
differences. It is found that cold galaxies are in regions with less 
"strangeness" than hot galaxies. 



If the model o f IRAS galaxies by 
iRowan-Robinson & Crawford (11989) is reasonable, the 
colour defined in the paper is closely associated with the 
SFR of galaxy. Effects of environment on star formation 
activity are very complex under competing mechanisms. In 
rich environments, spirals exhibit significant gas deficiencies 
resulting from ram pressure from intracluster mediums, 
together with possible tidal stripping from interactions with 
other galaxies and cluster potential; in this way, star formation 
rate is effectively suppressed. On the other hand, there are 
external environmental influences that can have much stronger 
effects on boosting the star formation rate. Studies of Ha and 
FIR emission of interacting and merging galaxies have shown 
strong excess star formation (see review of KennicuttL 119981) . 
Numerical simulations also pointed out pair interaction, and 
merging can produce tidal gas inflows during orbit decay; 
the inflow will drastically boost the star formation (e 



ally 



iBarnes & HernauistL[l995lTisseraLEo00llTissera et all 12002 

As displayed in Fig.0 in all subsamples, the possibility of 
finding an isolated galaxy is much higher than having galaxies 
with neighbours, while the "hot" subsample has a larger portion 
of galaxies with no companion than the "cold" subsample. The 
poor possibility of finding a galaxy with prominent present star 
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0.05 OJO 0J5 

Normalized counts of neighbours 

Fig. 7. Distribution of the number of neighbours of galaxies in 
different colour subsamples; the inset plot is the log-log version 
for a better demonstration, and the vertical arrow in the subplot 
marks the average normalised neighbour count for a randomly 
selected galaxy belonging to the main sample. The radius of the 
cell is r = 3.77/i"'Mpc. Counts of neighbours are normalised 
by the number of all galaxies, which is 1 1463, here. The num- 
ber of cells is counted within a bin of width 0.044, then it is 
normalised for each colour subsample by the total number of 
galaxies of that subsample. Error bars are estimated over 20 
bootstrap subsamples. 



formation activities in a rich environment can be interpreted 
as follows: in crowded regions, most galaxies experienced an 
active star formation stage on a very short time scale at high 
redshift, and then were quenched afterwards. The galaxies that 
were stiU there with considerable present SFR are likely the 
small number of survivors. This conclusion is similar to the 
analys is of the 2dFGRS and SDSS surveys by Balogh et alJ 
(120041) . Actually, observations of Spitzer have shown that the 
bulk of star formation in massive galaxies occurs at early cos- 
mic epochs and is largely complete by z ~ 1.5, while at z < 1 
low er-mass galaxies dominate the overall cosmic mass assem- 
bly (Papovich, 2006). It is highly possible that the IRAS galax- 
ies in rich regions a re of small mass, w hich is consistent with 
the discovery of IMartinez et al.l (120021) . Another aspect is that 
the field IRAS galaxies are probably systems of slow evolution, 
which may be tested by studies on the SFRs of galaxies with a 
"loose" friend and galaxies without a neighbour over a period 
of time. 

Nevertheless, in Fig. with the increasing number of 
neighbours, the probability of a hot galaxy living with such 
richness is smaller than that of a cold galaxy. This picture 
shows the downside effects of environment on the SFRs of 
galaxies, which is consistent with recent discoveries that high 
SFR galaxies are inhibited in rich environments in the 2d FGRS 
and SDSS surveys (Martinez et al., 2002: Dominguez et al.L 
l2002HBalogh et al.Ll2004 . 

In the high richness regime, the distribution of hot and 
cold galaxies are asymptotically in agreement with each other 



within error bars. There must also be some mechanism that ef- 
fectively enhances the galaxy's star formation rate, partly in 
compensation for the general suppression when a region be- 
comes very crowded, since otherwise the distribution curve of 
the "hot" s ubsample should be a lways lower than the "cold" 
subsample. iKrongold et al.l ( 120021) found that interacting IRAS 
galaxies have lower S loo/'^eo than isolated galaxies. It has also 
been found that in very dense regions, there is induced star for- 
mation in galaxy pairs at a very small separation, due to in- 
teractions, and that a galaxy with a high star formation rate 
is likely to trigger more star formation in its close compan- 
ion (e .g., Cartere^alj[_2001; Tissera et al., 2002. iLewis et all 
l2002t iSol Alonso et all 12004. .2006.) . However, in this regime, 
our distribution curves have large uncertainties as there are not 
many cells within the bin for counting since IRAS galaxies 
avoid rich clusters. We shall be conservative about this claim, 
which needs confirmation from future works with more com- 
plete and denser samples than the PSCz. 
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